Time-dependent Generalized SIC-OEP formalism and Generalized SIC-Slater 

approximation. 
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We propose a simplification of the full "2 sets" Time dependent Self Interaction Correction (TD- 
SIC) method, applying the Optimized Effective Potential (OEP) method. The new resulting 
scheme is called time-dependent "Generalized SIC-OEP". A straightforward approximation, using 
the spatial localization of one set of orbitals, leads to the "Generalized SIC-Slater" formalism. We 
show that it represents a great improvement compared to the traditional SIC-Slater/KLI formalisms. 

PACS numbers: 71.15.Mb,31.15.Ew,31.70.Hq,33.80.Eh 



Density Functional Theory (DFT) [l|, 0] has become a 
standard theoretical tool for the description of electronic 
properties in a variety of physical and chemical systems, 
especially for sizable numbers of electrons. The extension 
of DFT to Time-Dependent (TD) situations ^ IS a more 
recent achievement which also motivates numerous for- 
mal and practical investigations TDDFT contitutes 
to date one of the few, well founded theories, for describ- 
ing dynamical scenarios in complex systems. From the 
point of view of applications TDDFT requires simple ap- 
proximations, the simplest one being the Adiabatic Lo- 
cal Density Approximation (ALDA), which proved very 
useful in calculations of structure and low-amplitude ex- 
citations (optical response, direct one-photon processes) 
But self-interaction error plagues ALDA, which in 
particular spoils the dynamical description of excitations 
involving ionization, especially when close to threshold. 
A correct treatment requires a self-interaction correction 
(SIC), such the one proposed in This SIC has been 
used since then at various levels of refinement for struc- 
ture calculations in atomic, molecular, cluster and solid 
state physics, see e.g. [1, 0, H 0- But the SIC scheme 
leads to an orbital dependent mean field, which causes 
several formal and technical difficulties [10|, becoming 
all the more cumbersome in the TD case. Applications 
of SIC in TD situations arc thus mostly done in approxi- 
mate manner, linearized [ll| , using averaged-density SIC 
or relying on various approximations of TD Opti- 
mized Efl:ective Potential (TDOEP) scheme [H] (the lat- 
ter allows to deal with a common local potential). The 
most popular approximate OEP scheme is the Kriegcr- 
Li-Iafratc (KLI) approach 14, 1^, 1^ and, in a further 
step of simplification, the Slater approximation [13] ■ But 
both suffer from inconsistencies, such as violation of zero 
force theorem and energy conservation [l8l | , even if some 
applications may be more forgiving [Toj . 

We have recently generalized the considerations of sta- 
tionary SIC to the time dependent case [lOl and 
showed that the time propagation of the "exact" SIC 
scheme can be performed thanks to a double basis set for- 



mulation (both sets being connected by a unitary trans- 
formation building the same total density). The result- 
ing scheme satisfies all key formal properties. But the 
remaining potential becomes non-local, which substan- 
tially slows down numerical calculations. In the station- 
ary case, a SIC-OEP scheme relying on this two sets 
technique can be performed leading to "Generalized SIC- 
OEP" either in full formulation ^ or in the "Gen- 
eralized SIC-Slater" (GSlat) approximation [111 which 
emerges naturally because of the localization of one set of 
orbitals [23 . GSlat is numerically efficient and represents 
a great formal and practical improvement over traditional 
stationary SIC-Slater/KLI [H, H^l- It remains to extend 
Generalized SIC-OEP to the time domain and consider 
simple efficient approximations thereof, such as possibly 
the time-dependent GSlat (TDGSlat). It is the aim of 
this paper to address both aspects. 

The starting point is the SIC energy functional for elec- 
trons (all sums run over occupied states only ; we omit 
the (r, t) dependencies when it is not misleading and we 
note p = J2a IV-cP) 



^'sic = 



(1) 



XI I ^ I V'a ) + ^^ext [p] + i^ALDA [P] - X! ^ALDA [ I V'a P ] • 



The first term is the non-interacting kinetic energy, 
Ecxt [p] = / drpvcxt where Vcxt stands for the interac- 
tion with the ionic background and any other local pos- 
sibly time-dependent external field, and EaldaIp] is a 
standard ALDA energy-density functional The last term 
corresponds to the SIC. We mention in passing that the 
SIC, and with it all our following development, does also 
apply to more general funtionals as, e.g., the Generalized 
Gradient Approximation (GGA) The TDSIC equa- 
tions result from a variational principle applied to the ac- 
tion integral (here we use the standard action functional. 
Its causality and boundary conditions problems are dis- 
cussed in [3, [ill)- Furthermore, we explicitly include 
the orthonormality constraint with Lagrange multipliers. 
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This yields 



-5 Usic- / dtY,i^M>^^fA =0 



Asic = / ' dt(Esic - ^(V'a|ift9t|^a)) , (2) 

leading to one-body equations that we recast as [20t 

ihdtl^i) = hsiclVi) (3) 
Vi : = (V'/sIC/alda [IV'/aH - C/alda [l^al'] IV^cM 

where C/alda [Pq] = ('5-Ealda[p]/<5p)|p=p„ and 

/iSIC = /JALDA - X^^ALDA [IV'aP] |V'a)W'a| (5) 



^ALDA = 7; 1- Woxt + C/ALDA [p] ■ 

Note that /isic is applied to the diagonal orbitals \Lpi) 
in Eq. ([3]), but is calculated with the localized orbitals 
I ■0a) in Eq. ([5]). The |(/3i)'s are related to the |?/'q)'s by a 
unitary transformation 



(6) 



The wave functions (pi and the transformation coefficients 
Uia are varied independently where variation with respect 
to (fi yields Eq. ^ and variation of the Uia then leads 
to the symmetry condition ^ [2^, [23| . The symmetry 
condition guarantees conservation of orthonormality and 
enforces localization of the states IV'a) [l^l- It thus plays 
a key role in this formalism. Although Eq. ([3]) has proven 
to be tractable numerically [l^, [27|, the corresponding 
Hamiltonian is non-local, see Eq. ([S]), which implies a 
very high numerical cost. 

To develop a simplified scheme while maintaining local- 
ity, we apply this double set formulation to the TDOEP 
scheme. We start from the action integral The pre- 
vious considerations on the "exact" SIC energy showed 
that it should be expressed with the localized if) a, which 
satisfy the symmetry condition ([4]). We stationarize the 
action, imposing that the diagonal ipi (linked by a uni- 
tary transformation to the ipa) satisfy a Kohn-Sham like 
equation with a local and state-independent potential Vb 
(although OEP orbitals are not exactly the same as the 
SIC ones, we employ the same notations to keep the pre- 
sentation compact) 



^ALDA - Vo(r, t) (pi{T, t) = ihdt^Piiy, t) 



(7) 



The orthonormality constraint is implicitely contained 
in (l7|). The result will thus be considered as a local 
approximation of the SIC Hamiltonian The opti- 

mized effective potential Vq{v) is found by the variation 



SA^ic/ 5VQ{T,t) ~ which employs variations with re- 
spect to the Lpi through the chain rule for functional 
derivatives. The independent variation of the coefficients 
in ths transformation ([S]) remains as before. We ob- 
tain [III 

t At' f dr'(K,(r',i')-<(r',0) 
J "'-00 J 

xK,{r, t; r', t')ip* (r', t')ip,ir, t) + c.c. = (8) 



+ 00 



K,{v,t;v\t')^-i J2 ^*jir,t)^j{r',t')9it^t') (9) 
1 6 



Vi{r,t) 



ipi{r,t) Sip*ir,t) 



LDA[|'0a|^](^') 
J —00 



= ^TTtE^^WC^alda [M (r,0V'a(r,t)(10) 

to be fulfilled together with the symmetry condition (|4]) 
which again results from variation of the Uia ■ This is the 
TD "Generalized SIC-OEP" formahsm. The new feature 
in this double-set TDOEP appears in the Vi in (fTO|) which 
now employs the localized ^pa, in accordance with exact 
TDSIC. The previous applications of TDOEP to SIC as 
found in [l^ used the action ^ written in terms of 
the diagonal ipi, which leads to well known pathologies, 
as the incapability of the traditional SIC-Slater/KLI ap- 
proximations to reproduce the spatial localization. Eq. 
([8]) can be rewritten equivalently as 

Vo = ^e{Vs + Vk + Vc} - ^m{VTDi + Vtd2} 
where Vg, Vk, Vc are defined as 



p 



^ 1^ V.(p,V|(^,|^) 



(11) 

(12) 
(13) 



V^{v,t) = [^f l'dr'{Vo{r',t')-v:{r',t')) 



x^*(r',t')A',(r,t;r',i'). (14) 



To those potentials, which also appear in the stationary 
case, one has to add purely time-dependent contributions 



i 



dt'{p,{t')\v,{t')\p,{t')) 

(15) 



where 3i = ^^{'Pi^fi — fi'^'Pi) is the current den- 
sity. Some straightforward manipulation with the uni- 
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tary transformation Q allow to rewrite 



(Fo(r', i') - C/ALDAil^Z-aPKr', r', t')- (16) 



Gen. Slater 
SIC 




Note that the "Generalized SIC-KLI" approximation Vk 
has not the form that might be intuitively expected 21 1 . 
Moreover, as {tfi\vi\Lpi) = {ifi\vi\ipi)*, we have, with (fT5|) . 
9m{^TDi} = 0. 

It turns out that the fuU TD "Generalized SIC-OEP" 
scheme is very costly numerically, even more than the 
"exact" TDSIC formalism because time integrals (mem- 
ory effects) appear explicitly in the definition of Vq. Thus 
we propose a strong simplification which emerges natu- 
rally from the fact that the remain spatially very lo- 
calized d, H^l- This means that, at all time, at a given r, 
one single ijja mostly dominates the other wave functions 
ipfi^a and thus 



J2 ^t/ALDA[|#|']V'« ~ C/ALDAilV-aH^a. (17) 

One finds that the Slater contribution Vs is dominating, 
although the other terms are generally non negligible. 
But in the TD "Generalized SIC-OEP", if we approxi- 
mate Vo ^ Vs = J2a ^^T^C^ALDAilV'aP]; the localization 
(i.e. Eq. (jl7p ) amounts to have almost vanishing Vk and 
Pi, see Eq. p^ . and thus almost vanishing Vq and Vtd2, 
see Eqs. ^ and Thus, the TDOEP result natu- 

rally reduces to 



Vo 



E 



-C^alda[|V' a I 



(18) 



Eq. ([7]) generates the set ipi of diagonal states, which 
can be interpreted (to first order) as single electron or- 
bitals, while the unitary transformation ^ serves to ac- 
commodate the symmetry condition ^ which, in turn, 
defines the localized orbitals ipa entering the potential 
Vq as given in Eq. Note that this equation has the 

form of a Slater approximation [2§| but is constructed 
from the localized ipa and applied to the ipi. This is 
the time-dependent extension of the "Generalized SIC- 
Slater" scheme (TDGSlat). We showed in ^ that its 
stationary counterpart solves many problems encoun- 
tered with traditional SIC-Slater/KLI methods. 

We first check the performance of TDGSlat as com- 
pared to full TDSIC in Fig. [1] Times are shown relative 
to the time needed to propagate the cheapest solution 



FIG. 1: Gained time in TDGSlat compared to TD- 
SIC (in unit of time for ALDA propagation in Hn: 
t{Hn^^'^^'^°'^^) /tCRi^^^^^), and logarithmic scale); calculations 
on H chains as indicated. 



(ALDA) for Hydrogen chains on a given physical time 
interval (mind the logarithmic scale). The gain as com- 
pared to TDSIC is dramatic (typically of an order of 
magnitude) for a cost an order of magnitude larger than 
ALDA, for Hg. 
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FIG. 2: Total energy and dipole moment (along perturbation) 
in a quantum dot as a function of time in the linear domain. 
See text for details. 



It is well known that standard (one set) SIC-Slater 
and KLI approximations suffer from inconsistencies in 
the time domain. They can strongly violate energy con- 
servation and zero force theorem (ZFT) while nei- 
ther fuU TDSIC nor full "Generalized SIC-OEP" do 
suffer from these defects, but at the price of very heavy 
numerical cost. Strictly speaking, at TDGSlat level one 
has (for a non-explicitly TD external potential) 



dtEsic = — 3m / dr 
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FIG. 3: Same as Figure [2] but for a Nas cluster. 



ALDALIVqI 



dt^ifilpl^Pi) = / druextVp + ^Ke / dr 

i •' a •' 

x{[/ALDA[|^a|'] ~ ^ ^I/alDA [|# |'] jV'a A^- 

which in principle leads to energy and ZFT violation as 
well. However, the localization of the il)a, which implies 
the approximate relation p7|) . lets us hope that TDGS- 
lat produces "less" violated energy conservation and zero- 
force theorem, where "less" means that the conservation 
laws remain stable for a longer time span than in simple 
SIC-Slater/KLI. This alone could be progress for many 
practical purposes. It is thus a key issue to explore how 
TDGSlat practically performs in the time domain. This 
has been done in a variety of systems and we illustrate 
our results in Figures [2] and [3] on two typical cases: a 
quantum dot with 6 electrons, in the spirit of poj . and 
a small metallic cluster. Calculations have been done in 
full 3D using the same numerical methods as in [23| . In 
both examples we plot the time evolution of the total 
energy and the dipole moment after a small boost of the 
initial (ground state) electronic distributions which sim- 
ulates a very short laser pulse, and still allows to check 
energy conservation in time because the excitation field is 
switched off during propagation. We take care of consid- 
ering sufficiently moderate perturbations to remain in the 
linear domain (excitation energy a few percent of the typ- 
ical electronic level spacing) and we follow the dynamics 
over "long" times, typically over 100 eigenperiods of the 
system, as seen from the dipole oscillations. We restrict 
the calculations to TDGSlat and SIC-Slater. Full SIC 
calculations on such long times are prohibitively costly. 
We have checked on shorter times that they deliver per- 



fect energy and ZFT conservation. 

Both figures deliver the same message and we thus dis- 
cuss them together. While both fuU SIC and TDGSlat 
dynamics remain remarkably stable in time, standard 
SIC-Slater exhibits a sizeable violation of energy and 
ZFT (seen through the unstability of the dipolar mo- 
ment [18[). In the case of energy one observes a strong 
drift in time while for the dipole moment oscillations 
become extremely large (much larger than the original 
ones). When pursuing the TDGSlat over even longer 
times one observes a small energy violation, as expected 
on formal grounds. Still the effect is rather small and 
much delayed as compared to standard SIC-Slater. And 
the most important point is that the drift occurs on sufR- 
ciently long times so that relevant physics can be studied 
for shorter times. One can thus conclude that dynamics 
in the linear domain is controlled in TDGSlat. remaining 
very close to full TDSIC over long times. 

To summarize, we started from the two basis set formu- 
lation of TDSIC and apphed the TDOEP formalism to 
recover locality, resulting to the time-dependent "Gen- 
eralized SIC-OEP" formalism. As it is very costly nu- 
merically, we looked for a much less costly relevant ap- 
proximation, which naturally comes from the localized 
character of the V'q set and is called time-dependent 
"Generalized SIC-Slater" . By virtue of the double-set 
technique, it has a wider range of applicability than tra- 
ditional SIC-Slater/KLI approximations. In particular, 
we checked formally and numerically on various systems 
(organic, metal, quantum dot) that TDGSlat will satisfy 
much better energy and ZFT conservation. 

This work was supported, by Agence Nationale de 
la Recherche (ANR-06-BLAN-0319-02), the Deutsche 
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